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The  heat  dissipation  systems  which  have  liquids  as  heat  carriers  outperform  conventional  dissipation 
systems  at  TEGs  (thermoelectric  generators).  However,  new  elements  need  to  be  introduced  such  as 
pumps,  secondary  heat  exchangers  and  piping. 

A  predictive  computational  model  of  a  dissipation  system  involving  refrigerant  liquids  has  been 
implemented.  The  accuracy  of  the  model  is  93%  for  all  its  outputs:  the  total  thermal  resistance,  the 
hydraulic  losses  and  the  auxiliary  power  consumption.  The  validation  of  the  model  has  been  done  with  a 
prototype  mainly  composed  by  a  multi-channel  heat  exchanger,  a  fan-coil,  a  pump  and  several  sensors: 
temperature,  pressure  and  flow  meters. 

A  study  on  the  influence  of  the  water  and  the  air  mass  flow  over  the  total  thermal  resistance  has  been 
conducted.  The  total  resistance  dependence  on  the  air  mass  flow  shows  the  importance  of  including  the 
secondary  heat  exchanger  into  the  thermal  and  hydraulic  calculations.  The  smallest  resistance  does  not 
always  obtain  the  highest  net  power  generation,  the  high  demanding  power  of  the  auxiliary  equipment 
needed  to  obtain  this  resistance  influences  negatively  on  the  net  power  generation.  Among  the  exper¬ 
imental  points,  the  optimum  scenario  obtains  a  40%  additional  power  generation  with  respect  to  the 
smallest  resistance  point. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Compactness,  robustness,  reliability  and  lack  of  moving  parts 
are  the  main  advantages  of  TEGs  (thermoelectric  generators) 
(TEGs).  These  features  cause  less  maintenance  and  easiness  of 
control,  making  TEGs  a  better  solution  than  other  energy  conver¬ 
sion  systems  such  as  turbines  or  thermal  engines.  These  benefits 
did  not  get  unnoticed  by  the  aerospace  field;  thermoelectricity  has 
been  present  over  40  years  at  several  spacecraft  and  satellites,  with 
a  special  remark  on  Voyager  probes  capable  of  producing  150  W 
(-350  W/cm2)  with  7%  efficiency  [1],  The  latest  thermoelectric 
experience  in  the  outer  space  is  the  vehicle  space  mission  to  Mars 
[2]. 

On  the  earth,  several  researches  study  the  possible  application 
of  thermoelectricity  in  residual  energy  conversion.  The  origin  of  the 
waste  energy  can  be  very  wide:  vehicle  exhaust  gases,  electricity 
generation  or  conventional  industry  chimneys.  The  50%  of  the 
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power  consumption  of  automotive  electronics  can  be  generated 
recovering  part  of  the  waste  energy  being  thrown  to  the  atmo¬ 
sphere  [3],  a  10%  fuel  efficiency  increase  can  be  obtained  locating 
TEGs  on  the  exhaust  tailpipe  of  vehicles  [4], 

Despite  the  costless  nature  of  waste  heat,  the  reduced  efficiency 
that  thermoelectricity  presents  is  a  decisive  aspect  over  the  appli¬ 
cability  and  profitability  of  big  scale  systems.  Great  efforts  are  being 
made  to  improve  the  overall  efficiency  of  TEGs.  Two  are  the  main 
objectives:  the  optimization  of  materials,  in  order  to  improve  the 
figure  of  merit,  and  the  study  of  new  heat  exchangers  towards  the 
reduction  of  the  thermal  resistances  on  both  sides  of  the  TEMs 
(thermoelectric  modules).  The  importance  of  reducing  the  thermal 
resistance  between  the  heat  source  and  the  hot  side  of  the  TEM  and 
its  analogous  in  between  the  cold  side  and  the  ambient  was 
demonstrated  by  Astrain  et  al.  [5], 

A  high  heat  flux  per  area  is  what  makes  thermoelectricity  spe¬ 
cial  in  terms  of  refrigeration.  This  issue  occurs  in  other  applications 
such  as:  electronics,  laser  diodes  or  microchemical  reactors.  One  of 
the  solutions  focuses  on  heat  exchangers  provided  with  fluids  as 
heat  carriers  (involving  one  or  two  phases).  The  dissipator  can  be 
composed  by  small  channels  or  porous  media  [6],  For  example,  a 
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Nomenclature 

hydraulic  losses  kPa 

Ay 

absolute  measurement  error 

A  area  m2 

M 

fluid  dynamic  viscosity  kg/ms 

Cp  specific  heat  J /kgK 

D  pipe  diameter  m 

v 

fluid  kinematic  viscosity  m2/s 

/  friction  losses  coefficient 

Subscripts 

g  gravity  m/s2 

Gr  =  s^(Ts-Tm)L3  gj-^hof  number 

1 

TEM-to-water  heat  exchanger  entrance 

2 

TEM-to-water  heat  exchanger  exit 

h  convective  coefficient  K/m2W 

3 

water-to-ambient  heat  exchanger  entrance 

k  thermal  conductivity  W/mK 

4 

water-to-ambient  heat  exchanger  exit 

kiocai  local  losses  coefficient 

air 

air 

L  pipe  length  m 

amb 

ambient 

m  mass  flow  kg/s 

C 

thermoelectric  module  cold  side 

N  number  of  fins 

cond 

conduction 

Nu  =  nusselt  number 

Pr  “  k  Prandtl  number 

cont 

contact 

conv 

convection 

Q.  heat  power  exchanged/transmitted  W 

CP 

cold  plate 

R  thermal  resistance  K/W 

e 

exterior 

Re  =  reynolds  number 

exp 

experimental 

Ra=GrPr  rayleigh  number 

/ 

fan-coil 

T  temperature  K 

fin 

fin 

U  global  heat  transfer  coefficient  K/Wm2 

i 

interior 

v  fluid  velocity  m/s 

net 

net  power  generation 

W  generated/consumed  electric  power  W 

P 

pump 

sim 

computationally  simulated 

Greek  symbols 

sys 

piping  assembly 

p  fluid  density  kg/m3 

TEM 

thermoelectric  module 

e  material  roughness  m 

tot 

total 

77  efficiency 

water 

water 

thermosyphon  loop  cooling  system  obtains  double  refrigeration 
power  than  a  traditional  air  convection  system  [7], 

In  terms  of  thermoelectricity,  Zhou  et  al.  [8]  stated  that  liquid 
refrigerant  heat  exchangers  generate  higher  net  powers  than  con¬ 
ventional  finned  dissipators.  Nevertheless,  heat  exchangers  that 
include  liquids  need  pumps  to  circulate  the  refrigerant.  Each  heat 
exchanger  has  a  unique  flow  rate  that  obtains  maximum  net  power 
generation  [9],  At  low  mass  flows,  the  pumping  power  can  be 
negligible;  however,  the  opposite  occurs  when  the  mass  flow  is 
considerable  or  the  pressure  drop  is  significant.  Increasing  100 
times  the  pressure  drop  (from  0.5  to  50  kPa)  produces  an  increase 
of  3656  times  of  the  pumping  power  in  a  specifically  built  heat 
dissipation  system  for  a  TEG  [10], 

The  geometry  optimization  for  the  latter  heat  exchangers  fo¬ 
cuses  on:  number  and  geometry  of  the  channels,  flow  distribution 
or  the  inclusion  of  inserts.  The  reduction  of  the  diameter  of  the 
channels  leads  to  heat  transfer  enhancement  [  11  ],  parallel  channels 
obtain  higher  net  generation  that  serpentine  ducts  [12]  and  panel 
inserts  can  achieve  a  50%  net  gain  relative  to  the  absence  of  inserts 
[13],  Latest  researches  present  a  10%  enhancement  in  power  gen¬ 
eration  by  introducing  stirred  flows  into  the  heat  exchangers  [14], 
They  also  propose  the  use  of  high  thermal  conductivity  nanofluids 
as  a  feasible  optimization  field  [11], 

In  the  previous  mentioned  studies  there  is  no  treatment  of  the 
refrigerant;  the  costs  of  reusing  the  refrigerant  are  not  included. 
Water  can  be  used  as  refrigerant,  but  normally  there  are  no  big 
reservoirs  near  the  generators  where  the  water  can  be  taken  from. 

In  the  present  work  each  element  involved  in  the  refrigeration 
system  has  been  taken  into  account.  The  dissipator  located  on  the 
cold  side,  a  secondary  heat  exchanger  needed  to  reduce  the  tem¬ 
perature  of  the  refrigerant,  the  pump  and  the  piping  itself.  The 


computational  model  implemented  in  MATLAB  computes  the 
global  thermal  resistance,  the  hydraulic  losses  and  the  auxiliary 
equipment  consumption.  A  specifically  built  prototype  validates 
the  model  using  20  operating  points  that  present  different  values  of 
water  and  air  mass  flows.  This  global  computational  model  is  a 
powerful  optimization  tool  for  dissipation  systems.  Including  its 
outputs  in  the  TEG  computational  model  presented  in  Ref.  [5]  the 
net  generated  power  can  be  obtained. 

The  main  objective  of  this  study  is  to  create  and  validate  the  tool 
which  obtains  the  thermal  resistance  and  the  auxiliary  power 
consumption.  A  study  of  the  influence  of  water  and  air  mass  flows 
on  the  thermal  resistances  of  the  heat  dissipation  system  has  been 
included  as  well  as  a  sample  of  the  optimization  potential  that  the 
methodology  presented  in  this  paper  has. 

2.  Water-to-ambient  heat  dissipation  system  for  a  TEG  system 

A  water-to-ambient  heat  dissipation  system  for  a  TEG  dissipates 
the  heat  coming  from  the  TEM  to  the  heat  sink.  TEMs,  due  to  the 
Seebeck  effect,  convert  part  of  the  heat  received  on  their  hot  side 
into  electric  power.  The  rest  of  the  heat  is  emitted  through  the  cold 
side  to  the  heat  sink,  normally  the  ambient. 

In  this  work,  the  water  is  used  as  heat  carrier  due  to  its  high 
convection  coefficient.  Nevertheless,  the  dissipation  system  which 
involves  water  is  not  as  simple  as  in  Fig.  1 ;  the  water  needs  addi¬ 
tional  devices  such  as  pumps,  pipes  and  additional  heat  exchangers. 
Fig.  2  shows  every  element  of  the  dissipation  system  built. 

Below  the  cold  plate  the  TEMs  are  collocated.  Their  cold  sides 
are  in  direct  contact  with  the  multi-channels  heat  exchanger, 
named  cold  plate.  The  cold  plate  of  Fig.  2  dissipates  heat  form  the 
cold  side  of  the  TEMs  as  the  dissipator  of  Fig.  1  does.  However,  the 
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dissipation  system  of  Fig.  2  is  composed  by  two  heat  exchangers,  a 
pump,  pipes  and  valves  among  others.  Normally  TEGs  are  not  close 
to  rivers  or  water  reservoirs,  so  the  water  needs  to  be  cooled  down 
via  a  secondary  heat  exchanger. 

The  behavior  of  TEG  systems  is  simulated  by  Astrain  et  al.  [5] 
through  a  computational  model.  The  thermoelectric  net  power 
generation  is  the  output  of  this  model  while  the  thermal  resistances 
and  geometries  of  the  elements  of  the  TEGs  are  the  inputs.  This 
model  presents  95%  accuracy.  The  computation  of  the  total  thermal 
resistance  that  the  dissipation  system  of  Fig.  2  has  and  obtaining 
the  auxiliary  equipment  consumption  leads  to  obtain  the  thermo¬ 
electric  net  power  generation.  The  optimization  can  be  performed 
with  the  interaction  of  the  two  computational  models  presented 
above. 


but  also  on  the  air  mass  flow.  The  secondary  heat  exchanger,  which 
is  in  charge  of  reducing  the  water  temperature,  has  fans  that  help 
air  to  circulate  over  the  fins.  This  fact  is  showing  the  importance  of 
including  each  element  of  the  dissipation  system  into  the  calcula¬ 
tion  of  the  total  thermal  resistance. 

The  computational  model  is  composed  by  two  stages,  first,  the 
temperature  of  the  entries  and  exists  of  the  two  heat  exchangers 
are  obtained  to  secondly  solve  the  hydraulic  losses.  The  hydraulic 
losses  are  dependent  on  the  temperature  of  the  fluid;  this  is  why 
the  temperature  is  solved  in  a  first  instance. 

The  heat  flux  that  the  water  is  gaining  at  the  cold  plate  is  the 
starting  point  of  the  computational  model.  This  heat  exchanger  is 
located  on  the  cold  side  of  a  TEG,  thus  this  heat  flux  is  the  heat  that 
the  TEMs  are  evacuating  through  their  cold  side.  The  heat  gained  by 
the  water  at  the  cold  plate  is  the  heat  that  the  rest  of  the  elements 
included  in  the  dissipation  system  (fan-coil  and  pipes)  need  to 
dissipate  in  the  steady  state.  A  fixed  water  mass  flow  is  used  at  each 
calculation. 

Rtot  =  Tc~/amb  (l) 

Eq.  (1)  obtains  the  total  thermal  resistance  of  the  dissipation 
system,  but  Tc  is  unknown.  The  thermal  resistance  of  the  cold  plate 
is  the  term  that  connects  Tc  with  the  temperatures  of  the  key  points 
of  the  system,  Eq.  (2).  Each  temperature  cannot  be  obtained  sepa¬ 
rately,  because  it  depends  on  the  rest  of  them,  so  an  equation 
system  needs  to  be  solved. 

(Tc-h)-(Tc-T2)/ Infect}) 

Rrn  =  - B - V - 7  (2) 


The  temperatures  are  needed  to  calculate  the  thermal  re¬ 
sistances  of  each  element.  This  calculation  is  done  by  the  global 
heat  transfer  coefficient  presented  in  Eq.  (3). 


3.  Computational  model 

The  developed  computational  model  solves  the  thermal  re¬ 
sistances  as  well  as  the  hydraulic  losses  of  the  dissipation  system.  A 
water  mass  flow  range  is  simulated  obtaining  the  dependence  of 
the  thermal  resistances  and  hydraulic  losses  on  the  water  flow.  The 
total  thermal  resistance  does  not  only  depend  on  the  water  flow, 


(3) 


The  three  adding  terms  of  the  denominator  correspond  to  the 
internal  convective,  the  conduction  and  the  external  convective 
resistances.  The  thermal  resistance  of  each  component  is  obtained 
via  the  global  heat  transfer  coefficient,  Eq.  (3),  but  the  cold  plate. 


Fig.  2.  P&I  of  the  dissipation  system. 


P.  Aranguren  et  al.  /  Energy  74  (2014)  346—3 


349 


The  internal  convective  coefficients  are  given  by  Eqs.  (4)— (6). 
Laminar  flows  (Edwards)  are  represented  by  Eqs.  (4)  and  (5)  taking 
into  account  whether  the  flux  is  developed  or  not.  Eq.  (6)  dictates 
turbulent  flows  (Dittus  Boelter)  [15],  The  fluid  is  being  refrigerated, 
so  that  the  Dittus  Boelter  term  “n”  equals  0.3.  The  Prandt  numbers 
for  water  are  all  inside  the  range  in  which  the  equation  is  applied 
(0.7  <  Pr  <  160)  and  at  each  iteration  the  temperature  difference 
between  the  surface  and  the  liquid  is  checked  to  be  less  than  6  °C. 

Niij  =  3.66;  Developed  region ,  Re<  104  (4) 


Nui=  3.66 

0.0668  (%)  ReDPr 

H - — - jjj;  Undeveloped  region,  Re<104 

1  +0.04  (^(Y)  ReDPr^) 

(5) 


NUi  =  0.023Re£8Pr°  3,  Re  >  104  (6) 

Every  fluid  property  is  dependent  on  the  temperature  of  the 
fluid. 

The  pipes  have  circular  cross-sectional  area,  so  that  Eq.  (7) 
stands  for  the  conduction  resistance.  The  fan-coil  is  a  circular 
pipe-continuous  fin  heat  exchanger  where  water  flows  in  both  di¬ 
rections  and  it  is  cooled  down  by  the  air.  The  cross-sectional  area  of 
the  pipes  is  also  circular,  thus  Eq.  (7)  is  used. 


ln(§r) 

cond=2^mr 


(7) 


To  calculate  the  external  convection  two  scenarios  have  been 
considered,  natural  convection  (Churchill  and  Chu),  Eq.  (8),  and 
forced  convection  (Churchill  and  Bernstein),  Eq.  (9).  The  velocity  of 
the  surrounding  air  is  included  into  the  computational  model. 


( 


Natural  convection  (8) 


Nue  =  0.3 

0.62ReiptJ 

+HYy 


(_Reo_\ 
V 28200 ) 


Forced  convection 


As  the  fan-coil  is  a  circular  tube-continuous  fin  heat  exchanger, 
the  global  heat  transfer  coefficient  does  not  correspond  with  Eq. 
(3).  The  fin  efficiencies  need  to  be  included  into  Eq.  (3). 


cold  side  temperature  of  the  TEM  and  the  temperature  of  the  water. 
Three  are  the  terms  that  compose  the  thermal  resistance,  the 
contact  between  the  TEM  and  the  cold  plate,  the  conduction  and 
the  convection  resistances. 

RCP  =  Rcont  +  Rcond  +  Rconv  (11) 


Rcont  =  0.005  K/W 


(12) 


Rcond  =  0.0011  K/W  (13) 

The  values  of  the  contact  and  conduction  thermal  resistances  of 
the  cold  plate  (Eqs.  (12)  and  (13))  correspond  to  the  experimental 
values  that  the  heat  exchanger  located  on  the  cold  side  of  the  TEG 
has.  The  Nusselt  conventional  expressions  are  not  valid  to  calculate 
the  convective  term.  The  mini  geometry  of  the  cold  plate  and  its 
water  mass  flow,  laminar  along  the  entire  flow  range,  create  the 
necessity  of  modifying  the  conventional  expressions.  The  Nusselt 
number  correlations  are  not  capable  of  predicting  the  heat  transfer 
phenomena  for  laminar  flows  [16],  In  this  work  an  expression  based 
on  Sieder  and  Tate  equation  has  been  modified  to  obtain  accurate 
results. 

The  heat  fluxes  that  each  element  of  the  system  dissipates  are 
calculated  through  Eqs.  (14)— (17). 


UfAf((T3  -  Tamb)  -  (T4  -  Tamb)) 


(14) 


U23A2*m  -  ^mb)  ~  (T3  ~  Tamb)) 

infefe' 


(15) 


n  U41A41((T4-Tamb)-(T,  -Tamb)) 
541 ' 


(16) 


Q.sys  =  0.23  +  Q41  (^7) 

The  pipes  are  not  isolated,  so  the  heat  dissipated  through  the 
piping(Qsys)  is  present. 

An  iterative  process,  starting  with  the  supposition  of  7i,  T3  and 
T4  is  proposed  to  achieve  the  temperature  of  the  system.  T2  is  ob¬ 
tained  through  Eq.  (18),  the  heat  absorbed  by  the  water  at  the  cold 
plate  (Qc)  is  the  input  of  the  computational  model. 

(,8) 

Eq.  (18)  associates  the  entry  and  exit  temperature  of  the  cold 
plate  with  the  heat  that  the  fluid  is  gaining.  Eqs.  (14)— (16)  are  used 
to  build  similar  equations  for  each  element.  Thereby  a  linear  system 
of  three  equations  with  three  unknowns  (T],T3  and  T4)  emerges.  The 
iteration  process  continues  until  each  temperature  difference  be¬ 
tween  two  consecutive  iterations  drops  below  a  tolerance,  in  this 
case  10  6.  At  this  point  the  hydraulic  computation  starts. 

The  total  hydraulic  losses  are  the  sum  of  the  losses  of  each 
element  of  the  dissipation  system. 

APtot  =  APcp  +  APSys  +  A  Pf  (19) 


The  calculation  of  the  cold  plate  thermal  resistance  is  more 
complicated.  In  Eq.  (2),  the  thermal  resistance  stands  between  the 


Each  of  the  hydraulic  loss  is  divided  between  the  primary  loses, 
fl/D,  and  the  local  ones  kiocai- 
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The  friction  coefficient  (f)  depends  on  the  water  mass  flow.  The 
Reynolds  number  establishes  whether  Eq.  (21 )  or  Eq.  (22)  is  used. 


The  thermal  resistance  of  each  element  and  the  consumption  of 
the  auxiliary  equipment  are  the  outputs  of  the  computational 
model  developed  in  this  section. 

4.  Prototype 


(21) 


/=# 


-  1 .8  -  log  (  ( IgH  +  (  3^3: ) 


Re  >2300 


(22) 


The  water  mass  flow  and  the  cross  sectional  area  is  necessary  to 
obtain  the  water  velocity,  however,  the  velocity  in  each  channel  of 
the  cold  plate  is  unknown.  The  cold  plate  is  a  multi-channel  and 
local  the  losses  in  the  entry  and  exit  regions  of  the  channels  are  too 
significant  to  be  underestimated,  so  the  hypothesis  of  constant 
mass  flow  along  each  channel  stands  far  away  from  reality.  An 
iterative  non-linear  calculation  (friction  factor  depends  on  water 
velocity)  is  employed  to  obtain  the  hydraulic  losses  at  the  cold 
plate. 

Some  specific  values  for  the  secondary  losses  are  as  follows: 
kelbow— 0.5,  kentry  0.5,  kexit— 1  kbranchflow—  1  klmeflow  0.25  [15], 
Additional  local  losses  are  added  to  the  simulation  model  ac¬ 
counting  for  every  element  of  the  system  such  as:  valves,  sensors, 
expansion  tank  and  purgecock. 

The  water  mass  flow  increment  leads  to  a  new  calculation  of  the 
thermal  resistance  and  the  hydraulic  losses.  The  water  mass  flow 
upper  limit  indicates  the  end  of  the  simulation.  Fig.  3  shows  the 
computational  model  proceeding:  “AT  stands  for  the  water  mass 
flow  upper  limit,  while  the  lower  limit  must  also  be  chosen  by  the 
user. 


The  prototype  built  to  validate  the  computational  model  can  be 
seen  in  Fig.  4.  The  main  elements  that  are  present  are:  the  cold 
plate,  the  fan-coil  and  the  pump.  However,  other  elements  are 
needed  to  ensure  proper  operation  of  the  system  and  safeness:  a 
security  valve  to  prevent  excess  interior  pressure,  a  purgecock  and 
an  expansion  tank  to  store  possible  air  present  in  the  water  mass 
flow.  The  interior  volume  of  the  tank  is  2  L.  The  heat  introduced  into 
the  system  comes  from  a  hot  plate,  a  copper  plate  filled  with  10 
circular  electric  resistances.  The  maximum  power  that  the  hot  plate 
can  emit  is  2  kW.  The  hot  plate  dimensions  are  similar  to  the  cold 
plate  working  area,  195  x  140  mm2.  The  assembly  between  these 
two  elements  ensures  that  every  face  of  the  hot  plate  is  completely 
isolated  but  the  one  in  contact  with  the  cold  plate,  thus  all  the 
emitted  heat  goes  into  the  cold  plate. 

The  heat  power  introduced  into  the  dissipation  system  and  the 
heat  dissipated  through  the  secondary  heat  exchanger,  the  fan-coil, 
are  calculated  through  the  water  temperature  incremental,  Eq.  (23) 
and  Eq.  (24).  The  specific  heats  are  evaluated  at  the  temperature 
mean  of  each  heat  exchange. 

Qc  —  rhCp(T2  -  T])  (23) 

Q f  =  mCp(T3  -  T4)  (24) 

To  measure  temperature  and  pressure  at  the  entries  and  exits  of 
the  two  heat  exchangers,  four  temperature  and  pressure  sensors 
are  present.  An  additional  temperature  sensor  to  monitor  ambient 


Fig.  3.  Schematic  of  the  computational  model. 
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temperature  is  also  installed.  Two  extra  temperature  sensors  are 
added  to  measure  the  temperature  of  the  hot  plate.  All  the  tem¬ 
perature  sensors  are  thermocouples  of  type  K  while  the  pressure 
sensors  are  piezo-resistive.  An  electromagnetic  flow  meter  mea¬ 
sures  the  water  mass  flow  circulating  along  the  dissipation  system 
and  a  thermo-anemometer  measures  the  velocity  of  the  air  at  the 
fan-coil.  The  position  of  every  sensor  can  be  seen  more  clearly  in 
Fig.  2  and  their  resolution  and  accuracy  in  Table  1. 

The  pump  installed  is  a  continuous  pump  which  offers  the 
possibility  of  choosing  between  four  pumping  levels,  each  of  them 
defined  by  a  characteristic  curve.  The  pump  is  a  Jabsco  59520  series. 
The  maximum  water  mass  flow  that  this  pump  can  provide  is  21 1/ 
min  with  a  power  consumption  of  21  W. 

The  cold  plate  is  a  multi-channel  heat  exchanger.  It  is  formed  by 
a  total  of  26  channels,  each  of  them  with  a  length  of  148  mm  and  a 
6.2  mm  hydraulic  diameter.  The  length  of  each  manifold  is  176  mm 
and  their  diameter  is  9  mm.  The  cold  plate  is  provided  by  two 
entries  and  two  exits.  Its  external  dimension  is  230  x  190  mm2.  The 
external  area  accommodates  the  screws  in  charge  of  ensuring 
sealing  between  the  two  sides  of  the  cold  plate. 

The  fan-coil  is  a  circular  tube-continuous  fin  heat  exchanger.  24 
circular  tubes,  10  mm  diameter  together  with  226  continuous  fins 
form  the  kernel.  Three  fans,JAMICONJF1225S2H,  are  responsible  of 
forcing  the  air  between  the  fins.  The  maximum  electric  power 
consumption  of  the  three  fans  is  16.5  W.  The  forced  air  cross 
sectional  area  is  200  x  408  mm2. 

In  order  to  obtain  equal  external  conditions  at  every  experiment, 
the  prototype  is  located  in  a  climatic  chamber  set  to  22  °C. 

5.  Results  and  analysis 

5.1.  Experimental  results 

A  total  of  60  experiments  have  been  performed.  The  water  and 
air  mass  flows  are  the  variables  of  these  experiments.  Four  different 
working  levels  for  the  water  and  five  for  the  air  have  been  tested, 
reproducing  three  replicas  of  the  20  possible  experimental  sce¬ 
narios.  The  rest  of  the  parameters,  the  ambient  temperature  and 
the  heat  flux  introduced  intro  the  dissipation  system,  were  con¬ 
stant  all  along  the  experiments. 


Fig.  5  presents  the  working  levels  of  the  pump  and  fan-coil.  The 
power  consumption  tendencies  with  respect  to  the  water  and  air 
mass  flows  are  used  to  quantify  the  decrease  of  the  thermal  re¬ 
sistances  as  the  power  consumption  increases. 

Table  2  presents  the  experimental  temperatures  of  the  60  ex¬ 
periments  made,  three  replicas  of  the  20  working  points  studied. 
The  three  replicas  of  each  working  point  represent  the  deviation  of 
the  measurement  instruments,  a  maximum  of  0.6  °C.  This  deviation 
corroborates  the  accuracy  of  the  temperature  probes,  ±0.5  °C  as 
Table  1  shows.  An  additional  column, Tc-Tamb,  is  added.  This  col¬ 
umn  represents  the  total  thermal  resistance  of  the  dissipation 
system,  the  most  appreciable  thermal  resistance  since  it  is  the  value 
that  represents  the  whole  refrigeration  system  located  on  the  cold 
side  of  the  TEG.  The  absolute  experimental  error  of  this  indirect 
parameter  (Tc-Tamb)  is  calculated  through  Eq.  (25)  [17]  and  the 
accuracy  values  of  the  measurement  probes  stated  in  Table  1.  The 


Table  1 

Resolution  and  accuracy  of  the  sensors  employed. 

Sensor  Resolution  Accuracy 

Temperature  0.1  °C  ±0.5  “C 

Pressure  0.01  bar  ±1%  measured  value 

Water  flow  meter  0.01  1/min  ±0.2%  measured  value 

Air  velocity  0.01  m/s  ±0.01  m/s  +  3%  measured  value 


Fig.  5.  Working  levels  of  the  Fan-coil  and  the  pump. 


352 


P.  Aranguren  et  al.  /  Energy  74  (2014)  346—358 


Table  2 

Experimental  and  computationally  simulated  values. 


Working  point  7j(°C) _ 

Pump  Fan-coil  exp  sim  Error 


3  5 


27.3  27.05  0.25 

27.3  0.25 

27.2  0.15 

27.7  27.45  0.24 

27.5  0.04 

27.6  0.14 

28.1  27.88  0.22 

28.0  0.12 

28.3  0.42 

29.0  28.91  0.09 

28.9  -0.01 

29.3  0.39 

30.4  30.56  -0.17 

30.3  —0.27 

30.8  0.23 

27.1  27.14  -0.04 

27.1  -0.04 

27.1  -0.04 

27.5  27.40  0.09 

27.4  -0.01 

27.5  0.09 

28.0  27.79  0.21 

27.8  0.01 

28.1  0.31 

28.8  28.85  -0.05 

28.6  -0.25 

29.1  0.25 

30.3  30.62  -0.33 

30.0  -0.63 

30.6  -0.03 

27.0  27.14  -0.14 

27.0  -0.14 

27.1  -0.04 

27.4  27.46  -0.07 

27.3  -0.17 

27.4  -0.07 

27.8  27.97  -0.17 

27.7  -0.27 

27.8  -0.17 

28.6  28.92  -0.33 

28.5  -0.43 

28.8  -0.13 

30.0  30.38  -0.39 

29.9  -0.49 

30.5  0.11 

26.8  26.88  -0.08 

26.8  -0.08 

26.9  0.02 

27.1  26.98  0.11 

27.0  0.01 

27.2  0.21 

27.6  27.38  0.21 

27.5  0.11 

27.8  0.41 

28.5  28.72  -0.22 

28.3  -0.42 

28.6  -0.12 

29.8  30.12  -0.33 

29.7  -0.43 

30.2  0.07 


H°  C) 


exp 

29.1 

29.2 

29.1 
29.5 

29.4 

29.5 
30 
29.9 

30.2 

30.8 

30.7 

31.1 

32.2 

32.2 

32.6 
29.5 

29.5 

29.6 

29.9 

29.8 

29.9 

30.4 

30.2 

30.5 

31.1 

30.8 

31.4 

32.7 

32.4 

32.9 

29.9 
29.9 
30.0 

30.3 

30.1 

30.3 

30.8 

30.7 

30.8 

31.5 
313 

31.7 

32.9 

32.8 

33.3 

31.3 

30.9 

31.4 

31.4 

31.2 

31.6 
32 

31.7 
32 

32.8 
32.6 

32.9 
34 
34 

34.4 


28.91 

29.32 

29.78 

30.71 

32.40 

29.57 

29.80 

30.19 

31.11 

32.99 

30.04 

30.33 

30.97 

31.79 

33.25 

31.24 

31.28 

31.65 

33.02 

34.36 


0.18 

0.28 

0.18 

0.18 

0.18 

0.22 

0.12 

0.42 

0.08 

-0.02 

0.38 

-0.20 

-0.20 

0.20 

-0.08 

-0.08 

0.02 

0.09 

-0.01 

0.09 

0.21 

0.01 

0.31 

-0.02 

-0.32 

-0.29 

-0.59 

-0.09 

-0.14 

-0.14 

-0.04 

-0.03 

-0.23 

-0.03 

-0.17 

-0.27 

-0.17 

-0.29 

-0.49 

-0.09 

-0.36 

-0.46 

0.04 

0.05 

-0.35 

0.15 

0.11 

-0.09 

0.31 

0.35 

0.05 

0.35 

-0.22 

-0.42 

-0.12 

-0.36 

-0.36 

0.04 


r3(°c) 

exp  sim 


29.1  28.90 


29.5 

29.2 

29.4 
29.9 

29.8 

30.1 

30.8 
30.7 

31 

32.2 

32 

32.4 


29.4 

29.7 

29.6 


29.30 

29.76 

30.69 

32.37 

29.55 

29.78 


30.3 

30.1 

30.3 

31.1 

30.7 

31.3 

32.4 

32.2 

32.7 

29.7 
29.7 

29.9 

30.1 
30.0 

30.2 

30.6 

30.5 

30.6 

31.3 

31.2 

31.4 

32.7 

32.6 

33.1 
31 

30.7 

30.9 

31.3 
31.0 

31.2 

31.7 

31.5 

31.7 

32.5 
32.1 

32.6 

33.9 

33.8 
33.8 


30.16 

31.09 

32.96 

30.01 

30.30 

30.94 

31.76 

33.21 

31.20 

31.24 

32.96 

34.29 


0.19 

0.09 

0.09 

0.19 

-0.11 

0.09 

0.13 

0.03 

0.33 


0.30 

-0.18 

-0.38 

0.02 

-0.26 

-0.26 

-0.16 

-0.09 

-0.19 

0.01 

0.13 

-0.07 

0.13 

0.01 

-0.39 

0.21 

-0.56 

-0.76 

-0.26 

-0.32 

-0.32 

-0.12 

-0.21 

-031 


-0.34 

-0.46 

-0.56 

-0.36 

-0.52 

-0.62 

-0.12 

-0.20 

-0.50 

-030 

0.06 


0.10 

-0.10 

0.10 

-0.47 

-0.87 

-037 

-0.40 

-0.50 

-0.50 


ri(°c) 

exp  sim  Error 

27.2  27.05  0.15 

27.1  0.05 

27.1  0.05 

27.6  27.45  0.14 

27.4  -0.06 

27.5  0.04 

28.0  27.88  0.11 

27.9  0.01 

28.1  0.21 

28.9  28.91  -0.02 

28.8  -0.12 

29.1  0.18 

30.4  30.57  -0.17 

30.2  -037 

30.6  0.03 

27.0  27.14  -0.14 

27.0  -0.14 

27.1  -0.04 

27.4  27.40  -0.01 

27.3  -0.11 

27.4  -0.01 

27.9  27.79  0.11 

27.7  -0.09 

27.9  0.11 

28.8  28.85  -0.05 

28.5  -0.35 

29.0  0.15 

30.2  30.63  -0.43 

30.0  -0.63 

30.4  -0.23 

26.9  27.14  -0.25 

26.9  -0.25 

27.0  -0.15 

27.3  27.46  -0.17 

27.2  -0.27 

27.4  -0.07 

27.7  27.97  -0.28 

27.6  -0.38 

27.7  -0.28 

28.6  28.92  -0.33 

28.5  -0.43 

28.7  -0.23 

30.0  30.39  -039 

29.9  -0.49 

30.3  -0.09 

26.7  26.88  -0.19 

26.8  -0.09 

26.8  -0.09 

27.1  26.98  0.11 

27.1  0.11 

27.1  0.11 

27.5  27.38  0.11 

27.4  0.01 

27.7  0.31 

28.4  28.72  -033 

28.2  -0.53 

28.5  -0.23 

29.9  30.13  -0.24 

30  -0.14 

30  -0.14 


13.61 

14.08 

14.58 

15.16 

16.93 

14.79 

14.92 

15.28 

15.85 

17.93 

15.28 

15.54 

16.36 

16.95 

18.33 

17.11 

16.87 

17.23 

18.58 

19.77 


Error  (%) 


2.79 
1.74 

-3.90 

0.86 

0.16 

1.13 

1.13 

1.79 
2.53 
1.90 
2.84 

-3.54 

0.12 


-131 

0.73 

1.20 

0.54 

0.54 

2.05 

2.05 

3.35 

0.63 

4.23 


0.43 

0.75 

1.39 

133 

-0.26 

1.02 


-2.25 

-0.01 

-1.21 

-031 

0.11 

-0.44 

-0.71 

-1.57 


2.19 
3.58 

3.20 
2.38 
4.27 

-1.83 


-1.40 

0.38 

1.13 


maximum  relative  experimental  value  for  the  temperature  differ¬ 
ence  is  5.34%.  The  maximum  relative  experimental  error  represents 
the  measurement  deviation  obtained  when  the  temperature  dif¬ 
ference  is  calculated  experimentally. 


As  the  computational  model  is  used  to  obtain  the  total  thermal 
resistance  and  the  consumption  of  the  auxiliary  equipment,  the 


influence  of  the  pump  and  the  fan-coil  consumption  on  the  thermal 
resistances  is  studied  based  on  the  experimental  data.  Eqs.  (1 ),  (2) 
and  (26)  are  used  to  obtain  the  main  thermal  resistances  involved: 
the  total  resistance  and  the  resistances  of  the  cold  plate  and  the  fan- 
coil.  The  deviation  of  the  experimental  data  can  be  seen  in  Fig.  6, 
Fig.  7,  Fig.  8,  Fig.  9,  Fig.  10  and  Fig.  11.  These  figures  represent  the 
three  main  thermal  resistances  in  function  of  the  water  and  the  air 
mass  flows,  and  the  hydraulic  losses  as  a  function  of  the  water  mass 
flow.  The  three  replicas  of  each  working  point  represent  the 
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experimental  deviation.  Every  experimental  result  stays  inside  de 
±5%  deviation  range  due  to  the  measurement  equipment  as  can  be 
seen  from  Figs.  6  to  11. 


The  thermal  resistance  of  the  fan-coil  is  dependent  on  both  the 
water  and  the  air  mass  flows.  Figs.  6  and  7  indicate  stronger 
dependence  with  the  air  mass  flow  than  with  the  water.  The  global 
heat  transfer  coefficient  is  composed  by  internal  and  external 
convection  and  the  conduction.  In  this  case,  the  internal  convection 
is  water-forced  while  the  external  convection  is  forced  with  air.  The 
convective  coefficient  of  the  water  is  higher  than  the  air  one,  so 
that,  the  external  convection  is  defining  the  heat  transfer.  The 
improvement  of  the  external  convection  results  in  a  marked  ther¬ 
mal  resistance  reduction. 

The  resistance  of  the  cold  plate  only  depends  on  the  water  mass 
flow  as  presented  in  Fig.  8.  The  air  mass  flow  causes  changes  in  the 
temperature  of  the  water.  The  Nusselt  number  is  dependent  on  the 
Reynolds  and  Prandtl  numbers.  The  Reynolds  number  is  a  function 
of  the  water  mass  flow  and  the  fluid  properties,  including  the 
Prandtl  number,  dependent  on  the  temperature  of  the  water.  As  the 
temperature  changes  are  not  too  significant,  the  thermal  resistance 
of  the  cold  plate  is  independent  with  the  air  mass  flow.  The  vari¬ 
ability  between  points  is  due  to  the  experimental  deviation,  a  ±5% 
as  it  was  stated  at  the  beginning  of  the  “Results  and  analysis” 
section. 

Fig.  9  shows  the  dependence  of  the  total  thermal  resistance  on 
the  water  mass  flow  while  Fig.  10  presents  its  air  mass  flow 


dependency.  The  values  of  the  total  thermal  resistance,  Figs.  9  and 
10,  are  significantly  higher  than  the  ones  of  the  cold  plate.  Opti¬ 
mizing  the  dissipation  system  does  not  have  to  solely  focus  on  the 
thermal  cold  plate  optimization,  but  on  the  optimization  of  the 
system  as  a  whole.  The  total  thermal  resistance  strong  dependence 
on  the  air  mass  flow  displays  the  necessity  of  taking  into  account 
the  fan-coil  for  thermal  calculations. 

Changing  the  working  level  of  the  pump  from  1  to  2  produces  a 
51%  increase  in  the  water  mass  flow  (from  0.0239  kg/s  to  0.0361  kg/ 
s)  and  a  power  increase  of  the  87%.  The  cold  plate  thermal  resis¬ 
tance  reduction  is  12%  while  the  total  resistance  decreases  a  10%, 
being  the  fan-coil  set  to  level  1  (3.6  W).  On  the  other  hand,  rising 
the  pump  level  from  3  to  4,  the  reduction  of  the  total  resistance  is 
the  3%,  a  9%  decrease  on  the  cold  plate  resistance.  The  increase  of  a 
65%  in  the  power  consumption  of  the  pump  leads  to  an  increase  of 
the  26.5%  in  the  water  mass  flow.  Consequently,  the  effects  on  the 
total  thermal  resistance  of  increasing  the  working  level  at  low 
levels  are  higher  than  in  the  upper  levels. 

The  fan-coil  and  the  total  resistance  decrease  18%  and  7% 
respectively  if  the  level  of  the  fan-coil  is  increased  from  1  to  2  (the 
pump  working  level  is  set  to  4,  13.12  W).  The  consumption  in¬ 
creases  a  62.5%  and  the  air  mass  flow  a  43%.  If  the  fan-coil  is  set  to 
level  5,  coming  from  level  4,  the  consumption  increases  a  27%,  the 
water  mass  flow  a  12%  and  the  resistance  reduction  is  the  2%  and  3% 
for  the  fan-coil  and  total  resistances  respectively.  The  tendency  of 
the  total  thermal  resistance  is  quite  linear  with  the  fan-coil  power 
consumption. 

Fig.  11  presents  the  losses  characteristic  curve  of  the  system 
(P1-P4),  and  the  characteristic  curves  of  different  parts  of  the 
dissipation  system.  Each  sensor  position  can  be  seen  in  Fig.  2.  The 
intersection  of  the  characteristic  curve  of  the  pump  with  the 
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hydraulic  losses  curve  of  the  system  obtains  the  operating  point, 
the  water  mass  flow  that  circulates  along  the  dissipation  system. 
Using  a  high  power  demanding  pump,  a  high  water  mass  flow 
could  have  been  achieved  obtaining  a  very  low  thermal  resistance. 
In  this  case,  the  power  consumption  would  be  very  high  leading  to  a 
non-desirable  solution,  lower  net  thermoelectric  generation. 

The  hydraulic  losses  of  the  fan-coil  (P3-P4)  are  the  smallest 
ones,  but  the  only  ones  that  just  take  into  account  a  single  element 
of  the  system.  The  hydraulic  losses  between  point  1  and  2  (P1-P2) 
take  into  account  the  cold  plate  itself,  the  piping  needed  to  circulate 
the  water,  a  purgecock  and  a  temperature  sensor.  The  hydraulic 
losses  between  point  2  and  3  (P2-P3)  represent  the  losses  due  to 
the  additional  elements  of  the  system,  an  expansion  tank,  a  security 
valve,  piping,  elbows  and  temperature  and  mass  flow  sensors. 
These  results  present  the  magnitude  of  the  local  losses. 

As  Fig.  11  shows,  each  hydraulic  loss  is  essential  to  compute  the 
total  hydraulic  losses.  However,  a  final  design  would  not  need  to 
include  so  many  additional  elements,  so  that  the  total  hydraulic 
losses  would  be  reduced  drastically. 

5.2.  Computational  model  validation 

To  obtain  the  simulated  values  the  experimental  water  mass 
flow  is  introduced  into  the  computational  model.  The  total 
computational  time  needed  to  calculate  the  60  different  cases  is 
1.15  s  employing  a  3.10  GHz  CPU  (Central  Processing  Unit).  It  takes  a 
media  of  0.019  s  to  calculate  the  thermal  resistances  and  the  total 
hydraulic  losses  of  a  particular  case. 

The  temperature  tendencies  are  well  reproduced  by  the 
computational  model  as  Fig.  12  and  Fig.  13  show.  The  absolute  error 


Water  mass  flow  (kg/s) 


Fig.  12.  Experimental  temperature,  rlexp  and  /X.xp,  and  computationally  simulated 
temperatures,  Tlsim  and  T2s jm  of  the  cold  plate. 


of  the  temperatures  is  obtained  through  Eq.  (27).  This  error  is 
bounded  by  its  minimum,  -0.87  °C,  and  its  maximum,  0.42  °C,  as 
can  be  seen  in  Table  2.  In  the  case  of  the  temperature  difference  that 
represents  the  total  thermal  resistance,  the  relative  error  has  been 
calculated,  the  range  where  each  error  can  be  found  is  [-3.90  to 
4.27%].  A  very  small  temperature  drops  exist  between  T2  and  T3  and 
between  T4  and  Ti,  the  areas  in  between  the  cold-plate  and  the  fan- 
coil  as  Figs.  12  and  13  present.  Fig.  2  schematically  presents  the 
elements  included  in  each  area;  pipes,  measurement  probes  and 
safety  elements.  The  characteristics  of  the  prototype  determine  the 
heat  emitted  to  the  atmosphere,  material  and  length  of  the  pipes. 
Nevertheless,  in  an  installation  where  the  length  of  the  pipes  is 
considerable  and  the  material  of  the  pipes  is  a  better  heat 
conductor,  the  heat  emitted  to  the  atmosphere  by  these  areas  can 
achieve  important  values. 

Absolute  error  =  ValueeXp  -  ValueSjm  (27) 


Relative  error  = 


Valueexp  -  Valuesim 
Valueexp 


100 


(28) 


Figs.  14—17  show  the  computational  model  accuracy.  Figs.  14—16 
present  the  computational  simulated  values  versus  the  experi¬ 
mental  values  of  the  resistances  of  the  fan-coil,  the  cold  plate  and 
the  total  one,  respectively.  As  can  be  seen,  every  simulated  value 
stays  in  a  ±7%  error  range,  thus  the  accuracy  of  the  computational 
model  is  93%. 

The  total  resistance  is  meant  to  be  introduced  into  the  TEG 
computational  model  to  obtain  the  net  power  generation.  Fig.  17 
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presents  the  statistical  study  of  the  error;  the  standardized  bias  and 
kurtosis  are  -0.3116  and  -0.3188  respectively,  so  the  assumption  of 
normality  can  be  made.  As  Fig.  17  b)  shows,  the  standardized  errors 
stay  in  between  the  -3  and  +3  range;  just  3  errors  are  higher  than 
|2|.  The  mean  of  the  error  is  0.22. 

Fig.  18  presents  the  total  hydraulic  losses  predictability  of  the 
computational  model.  The  computational  model  maximum  error  is 
±7%.  An  accurate  prediction  of  the  hydraulic  losses  is  vital  for  pump 
design. 

The  rapidity  of  the  computational  model  and  its  versatility 
converts  it  into  an  optimization  tool  to  dissipation  design.  Changes 
in  geometry,  auxiliary  equipment  or  ambient  conditions  can  be 
introduced  and  evaluated  seeking  net  power  generation 
maximization. 

The  computational  model  also  obtains  accurate  results  (a 
maximum  error  of  ±7%)  if  the  pump  characteristic  curve  is  intro¬ 
duced,  Fig.  20  depicts  this  fact.  The  pump  curve  is  determinant, 
introducing  an  overestimated  curve  would  lead  into  a  bigger  water 
mass  flow  operating  point,  so  that  a  smaller  thermal  resistance  than 
the  measured  value. 

An  accurate  water  mass  flow  prediction  is  vital  for  resistances 
prediction.  Fig.  19  provides  the  simulated  values  versus  the 
experimental  data,  the  accuracy  is  93%.  The  computational  time  to 
obtain  these  results  grows.  The  model  has  to  calculate  in  each  case 
the  outputs  in  a  range  of  water  mass  flows.  The  bigger  the  number 
of  calculated  points  inside  the  flow  range  the  better  results  are 
obtained.  The  total  computational  time  for  the  60  cases  is  26.26  s. 
This  value  was  expectable  because  each  case  needs  to  compute  26 


water  mass  flow  points  in  order  to  obtain  the  operation  point;  the 
exact  water  mass  flow  that  circulates  along  the  system. 

The  statistical  study  of  the  error  is  presented  in  Fig.  21.  The 
standardized  bias  and  kurtosis  are  -0.3429  and  -1.4276  respec¬ 
tively  and  the  standardized  errors  are  all  included  in  the  range 
between  -2  and  2  with  the  exception  of  a  particular  error  with  a 
value  of  2.10.  The  statistical  study  on  the  errors  assures  that  the 
model  is  valid  with  a  maximum  error  of  ±7%  regardless  of  whether 
the  mass  flow  or  the  pump  curve  is  introduced. 

6.  Optimization  study 

The  computational  model  presented  and  validated  in  the  above 
sections  procures  the  total  thermal  resistance  and  the  auxiliary 
power  consumption  of  a  heat  dissipation  system  involving  water  as 
a  heat  carrier.  In  thermoelectricity,  the  thermal  resistances  on  both 
sides  of  the  TEMs  are  important  to  be  as  low  as  possible,  but  at  the 
same  time,  it  is  essential  to  minimize  the  power  consumption  of  the 
auxiliary  equipment,  therefore,  a  compromise  should  be  reached. 
The  validated  thermoelectric  generation  computational  model  [5] 
is  used  to  conduct  this  study. 


6.1.  Computational  optimization 


The  interaction  of  the  two  validated  computational  models  is 
used  to  search  the  optimum,  the  maximum  net  generation  working 
point.  An  expression  of  the  total  thermal  resistance  of  the  prototype 
thanks  to  the  computational  methodology  presented  in  this  work 
has  been  obtained.  The  expression  is  bounded  by  the  upper  and 
lower  experimental  limits  of  the  water  and  air  mass  flows 
0.02383  <  rhwater  <  0.05633  and  0.0906  <  mair  <  0.2432. 


Rtot(niwater,  mair)  =  0.044  -  0.5782fhWater  -  0.0015  ln(mair) 

■*  4.2272 m^ater  +  0.0016  ln(mair)2 
-  0.01m water  ln(mair) 

(29) 


Eq.  (29)  presents  the  dependency  of  the  total  thermal  resistance 
on  the  water  and  air  mass  flows.  This  equation  is  introduced  into 
the  thermoelectric  generation  computational  model  as  the  thermal 
resistance  of  the  dissipator  of  the  cold  side.  A  temperature  of  250  °C 
for  the  heat  source,  a  ambient  temperature  of  15  °C  and  a  thermal 
resistance  of  0.033  K/W  for  the  heat  exchanger  located  on  the  hot 
side  have  been  simulated. 

The  thermoelectric  generation  is  presented  in  an  equation,  Eq. 
(30),  as  a  function  of  the  water  and  air  mass  flows.  The 
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thermoelectric  generation  is  maxima  at  the  upper  bounds  of  the 
water  and  air  mass  flows.  Eq.  (31 )  addresses  the  net  generation;  the 
consumption  of  the  auxiliary  equipment  has  been  subtracted  from 
the  thermoelectric  generation.  In  this  case  the  optimum  net  power 
corresponds  to  0.0319  and  0.1315  kg/s  for  the  water  and  air  mass 
flows  respectively.  The  net  power  production  of  the  optimal  point 
equals  41.63  W  Fig.  22  presents  both  generations  highlighting  the 
importance  of  introducing  the  auxiliary  consumption  into  the 
calculation.  Whereas  the  lowest  thermal  resistance  generates  the 
highest  thermoelectric  power,  the  optimum  net  generation  power 
is  close  to  the  opposite  point,  the  smallest  auxiliary  consumption. 


Experimental  hydarulic  losses  [kPa] 


Fig.  18.  Total  simulated  hydraulic  losses  in  the  +7%  experimental  losses  range. 


WTEM(>riwater,fiiair)  =  (o.0208  -  0.1603mwater  -  0.0003  ln(mair) 
+ 1 .2587m^,ater  +  0.0005  ln(mair)2 
+  0.001 9mwater  ln(mair))  ~ 

(30) 

VVnet(mwater,  mair)  =  23.0931  +  560.302mwater 

+146.205mair  -  9218.5m^ater 
-581.183mair  +  207.173mWatermair  (31) 

6.2.  Optimal  experimental  point 

The  same  study  as  in  the  previous  section  has  been  made.  In  this 
case  the  resistance  values  of  the  dissipator  located  on  the  cold  side 
are  experimentally  obtained.  Four  water  mass  flows  included  in  the 
range  (0.02383-0.05633)  and  five  air  mass  flows  (0.0906-0.2432) 
have  been  chosen  to  obtain  the  thermal  resistances.  These  twenty 
experimental  points  correspond  with  the  experiments  made  to 
validate  the  computational  model  developed  in  this  work.  The 
same  values  of  ambient  temperature,  heat  source  temperature  and 
thermal  resistance  of  the  dissipator  of  the  hot  side  have  been  used, 
15  °C,  250  °C  and  0.033  K/W  respectively.  The  thermoelectric 
generation  as  well  as  the  net  generation  can  be  seen  in  Fig.  23.  The 
optimal  working  area  corresponds  to  the  two  lowest  water  and  air 


Experimental  thermal  resistance  [K/W] 


♦  Rtot 


Fig.  20.  Total  simulated  thermal  resistance  in  the  ±7%  experimental  resistance  range 
including  the  pump  characteristic  curves  as  inputs. 


Fig.  19.  Simulated 


;  flow  in  the  ±7%  experimental 
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mass  flows  (0.02383  <rnwater<  0.03622, 0.0906  <mair<  0.1420). 
The  net  generation  is  approximately  40  W  for  each  point  in  this 
area.  It  can  be  seen  that  the  computational  and  experimental 
optimal  values  agree. 

Fig.  23  presents  the  thermoelectric  and  net  power  generations. 
As  the  pump  and  the  fan-coil  consumption  increase  the  thermo¬ 
electric  power  generation  increases  also.  However,  the  net  gener¬ 
ation  does  not  have  this  tendency.  The  consumption  of  the  auxiliary 
equipment  is  not  negligible,  being  decisive  for  optimization 
seeking.  The  net  thermoelectric  generation  of  the  optimal  point  is 
40.73  W  while  the  smallest  thermal  resistance  working  point,  the 
point  in  which  the  power  consumption  of  the  auxiliary  equipment 
is  the  highest,  has  a  power  generation  of  28.99  W.  Smaller  mass 
flows  of  air  and  water  obtain  an  increase  of  the  40%  with  respect  to 
the  best  thermal  scenario,  the  smallest  thermal  resistance  obtained 
through  the  highest  air  and  water  mass  flows. 

7.  Conclusions 

A  computational  model,  which  represents  a  whole  dissipation 
system  for  a  TEG,  has  been  created.  The  main  elements  that  this 
model  has  are:  a  multi-channel  heat  exchanger  located  on  the  cold 
side  of  the  TEM,  a  secondary  heat  exchanger  in  charge  of  reducing 
the  temperature  of  the  water,  a  pump  to  impulse  the  water  and 
pipes  to  connect  the  elements. 


Fig.  22.  Thermoelectric  and  net  generated  power  computationally  calculated. 


A  prototype  has  been  built.  The  prototype,  as  the  computational 
model,  includes  each  element  needed  in  a  closed  dissipation  sys¬ 
tem,  apart  from  the  necessary  measurement  instrument  needed  to 
collect  the  data. 

It  has  been  proved  that  it  is  necessary  to  include  the  secondary 
heat  exchanger,  the  fan-coil,  and  the  piping  into  the  thermal  and 
hydraulic  calculations.  The  total  thermal  resistance  and  the  hy¬ 
draulic  losses  are  not  only  determined  by  the  cold  plate,  but  also  by 
the  rest  of  the  elements  of  the  system,  the  fan-coil,  the  piping  and 
the  additional  elements.  So  that,  every  element  of  the  system 
should  be  included  into  the  calculations  to  obtain  accurate  results. 

An  87%  pump  consumption  increase  leads  to  a  10%  total  thermal 
resistance  reduction  at  low  pump  levels.  A  65%  rise,  at  high  pump 
levels,  obtains  a  thermal  reduction  of  the  3%.  On  the  other  hand,  a 
62.5%  increase  in  the  fan-coil  consumption  obtains  a  7%  total 
resistance  decrease  at  any  working  level. 

The  computational  model  accuracy,  introducing  the  exact  water 
mass  flow  value,  or  with  the  pump  characteristic  curve,  is  93%.  The 
statistic  study  of  the  errors  verifies  the  model.  The  total  thermal 
resistance  and  the  auxiliary  equipment  consumption  introduced 
into  the  validated  TEG  computational  model  [5],  obtain  thermo¬ 
electric  net  power  generations,  the  maximization  objective.  Thus, 
an  optimization  tool  for  dissipation  systems  involving  refrigerants 
has  been  designed. 


WTEM 


Fig.  23.  Thermoelectric  and  net  generated  power  obtained  from  the 


experimental 
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The  optimal  configuration,  the  configuration  that  obtains  the 
highest  net  power  does  not  match  with  the  configuration  with  the 
smallest  total  thermal  resistance.  The  power  consumption  of  the 
auxiliary  equipment  needed  to  obtain  small  thermal  resistances 
consumes  a  big  part  of  the  generated  power  leading  to  undesirable 
situations.  The  optimum  experimental  point  obtains  an  increase  of 
the  40%  on  the  net  generated  power  with  respect  to  the  smallest 
resistance  experimental  point. 
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